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Abstract. The Gunn-Peterson effect, absorption of Lya photons by a 
homogeneous component of the intergalactic medium (IGM), potentially 
provides a test of Big Bang Nucleosynthesis (BBN). With a lower limit 
on the UV radiation field estimated from the contribution due to QSOs, 
a measurement of the Lya opacity of the intergalactic medium would 
permit the derivation of a lower bound to the baryonic density of the uni- 
verse. The effect, however, has continually eluded a convincing detection, 
both in H I and He II, despite extensive searches. Recent cosmological 
^Cg ■ hydrodynamical simulations of structure formation in the intergalactic 

<^> medium suggest an explanation for its absence. In a Cold Dark Matter 

■ dominated cosmology, the fragmentation of the baryons is nearly com- 

plete, leaving a negligible remnant to comprise a smoothly distributed 
component. The fragmentation extends even into regions that are under- 
. dense, where it gives rise to most of the optically thin H I systems and 

nearly all of the He II systems, both thin and saturated. The result is a 
Lya opacity from a smooth IGM that is suppressed by over two orders 
i of magnitude from the BBN value. 
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1. Introduction 



Soon after the identification of the first QSO, Gunn & Peterson (1965) recog- 
nized that its spectrum could be used to place a stringent limit on the density of 
neutral hydrogen in the intergalactic medium (IGM). Photons blueward of the 
frequency of Lya will be absorbed by neutral hydrogen en route to the observer 
in an expanding homogeneous IGM while photons redward will be transmitted 
(Field 1959). Hence, a step in the spectrum of a QSO is expected shortward 
of its Lya emission line. Gunn and Peterson found a step corresponding to an 
opacity of ~ 1/2, and a neutral comoving H I density of nni ~ 2 x 10 -12 cm -3 . 
This is only a tiny fraction of the density of hydrogen in galaxies today. Their 
inference was that the IGM must be highly ionized. The detection of absorp- 
tion by a homogeneous ionized component, the "Gunn-Peterson effect", would 
have important implications for cosmology. If the gas is photoionized, then the 
equation of photoionization equilibrium may be solved for the baryon density of 
the IGM, given an estimate of the ionizing radiation field. A lower bound on 
the radiation field may be placed by summing the contribution due to QSOs, 
which would then permit a lower bound to be set on the density of the IGM. 
This would be a direct test of Big Bang Nucleosynthesis (BBN). 
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Subsequently, however, more detailed spectra revealed a host of absorption 
lines shortward of Lya. These were identified by Lynds (1971) as Lya absorp- 
tion by intervening discrete absorbers, the Lya forest. The original opacity 
measured by Gunn and Peterson could be accounted for by the forest. Because 
the forest is clumped, it is not possible to solve for the baryon density of the 
IGM independently of a cloud model. Searches since then for absorption by a 
homogeneous component of the IGM (e.g., Steidel & Sargent 1987), led to sim- 
ilar results: all the absorption could be accounted for by the Lya forest. The 
homogeneous unclustered intergalactic medium expected from Big Bang nucle- 
osynthesis has continued to evade detection. Recent numerical hydro dynamical 
computations of structure formation in Cold Dark Matter (CDM) dominated 
cosmologies suggest an explanation. These simulations show that all the baryons 
undergo condensation into discrete systems, even when the density of the system 
is less than the cosmic mean. These systems give rise to absorption lines that 
are optically thin at line center, and form in regions a few megaparsecs across 
that are themselves underdense, "cosmic minivoids." 

2. Baryons Lost 

The Lya opacity of homogeneously distributed intergalactic hydrogen photoion- 
ized by an ambient metagalactic radiation field is 

r a ~ o.i8h^T^ - 75 (n D h 2 50 ) 2 (i + zf{\ + 2 ft) zr 1/2 r III 1 i _ 12 , (1) 

where flr> is the total baryonic density of the diffuse homogeneous component 
in terms of the closure density, T4 is the temperature of the neutral hydrogen 
in units of 10 K, Fhi,-i2 is the photoionization rate of the neutral hydrogen in 
units of 10 _12 s _1 , qo is the cosmological deceleration parameter, and /150 is the 
Hubble constant in units of 50kms _1 Mpc -1 . In the redshift range 2 < z < 3, 
the contribution of QSOs to the photoionizing background gives 1^1,-12 ~ L 
with an uncertainty of perhaps a factor of 2-3 (Meiksin & Madau 1993; Haardt 
& Madau 1996). By z = 4 this value may decline by a factor of 2, and by a 
factor of 20 by z = 5, although the QSO counts become increasingly uncertain 
at these high redshifts. The baryon density due to BBN is estimated to lie in 
the range 0.04 <; f^io <; 0.08 (2a, Copi et al. 1995), consistent (though just 
barely!) with recent D/H measurements in intergalactic gas clouds (Tytler et 
al. 1996; Rugers et al. 1996). The expected IGM opacity is thus r a ~ 0.15 at 
z = 3 and 0.8 at z = 4. These values greatly exceed current (2a) limits. Steidel 
& Sargent (1987) find r a < 0.08 at z « 3, while Giallongo et al. (1994) find 
T a < 0.08 at z w 4.3, although perhaps a more conservative upper limit is that 
of Jenkins & Ostriker (1991) of r a < 0.6 at z « 4. 

The conflict between the expected and measured upper limits may be allevi- 
ated by postulating radiation sources in addition to QSOs, like stars or decaying 
particles, or if a substantial fraction of high redshift QSOs were missed in optical 
surveys due to obscuration by dust in intervening absorbers (Fall & Pei 1993). 
An increase by a factor of up to 10 would be consistent with the proximity ef- 
fect (Bajtlik et al. 1988), although smaller values derived from the proximity 
effect may be favored (Espey 1993). Alternatively, Meiksin & Madau (1993) 
suggested that the conflict may indicate a very clumpy IGM. They argued that 
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a large fraction of the baryons may be contained in the Lya forest, given the 
observed numbers and size constraints of the absorbers. Their estimate was 
0.002 < < 0.05 for the density parameter of baryons in the forest. 

3. Numerical Simulations 

It has recently become possible to compute the growth of structure in the IGM 
using combined N-body/hydrodynamics numerical simulations. Unlike galaxy 
formation, the evolution of the IGM is a relatively "clean" problem since metal 
cooling is negligible and feedback effects from star formation are arguably unim- 
portant, or at most only affect the ionizing radiation field. Within the context 
of a given cosmology and primordial power spectrum, it is thus possible to treat 
the calculation as an initial value problem, without the need for introducing 
additional physical processes "by hand." 

Several groups have performed simulations of the IGM using a variety of 
techniques, and all have had remarkable success in reproducing the principal 
observed properties of the Lya forest (Cen et al. 1994; Zhang et al. 1995; 
Hernquist et al. 1996; Miralda-Escude et al. 1996). An unexpected finding 
of these calculations is that the absorbers arise predominantly in sheets and 
filaments that form coherent structures over scales of a few megaparsecs, with 
a typical thickness of a few hundred kiloparsecs. These calculations readily 
reproduce the observed near power-law column density distribution for the forest 
over the range 13 < logA^i < 17. 

In this talk, I present results from a simulation of the IGM in a Cold Dark 
Matter dominated cosmology (Zhang et al. 1996a,b). Standard CDM is able 
to match the clustering properties of material from galaxy cluster scales down, 
although some variant of the model is required to match to the fluctuations 
measured by COBE on larger scales. The simulations are performed in a flat 
cosmology with no cosmological constant, and /150 = 1. The normalization of 
the power-spectrum is given by as = 0.7. Since the neutral density of individual 
clouds scales like Qbh^o/Tui, the counts of objects depend on the magnitude 
of the radiation field. The results presented below are for the radiation field 
determined by Haardt &; Madau (1996). The hydrodynamics is solved on an 
Eulerian grid, with a comoving resolution of 75 kpc on a top grid of comoving 
size 9.6 Mpc. The behavior of the gas is resolved in a second subgrid, centered 
in the top grid and with a resolution 4 times greater. The Eulerian nature of the 
calculation has the important advantage over Lagrangian schemes of being able 
to resolve structures in underdense regions. These structures will be shown to 
be crucial for interpreting the amount of absorption by diffuse H I and especially 
He II gas in the IGM. Further details of the simulation are provided in Zhang 
(1996) and Zhang et al. (1996a, b). 

4. Baryons Regained 

Figure 1 shows the column density distribution of the clouds from the simulation. 
The clouds are found using a deblending algorithm for the lines, essential to avoid 
undercounting weak optically thin systems (Zhang et al. 1996a). The column 
densities are computed by fitting Voigt profiles to the lines. The agreement with 
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Figure 1. (Left panel) The distribution of HI column density at z = 3. 
The dotted line marks the division between optically thin and optically 
thick systems at line-center, for b = 15kms _1 . (Right panel) The mean 
baryon density, cloud temperature, peculiar velocity divergence, and 
cumulative baryon distribution as a function of HI column density. 



the observed distribution is amazingly good, as was found by the other groups 
for the optically thick systems, having H I column densities logAni > 13.3 for 
a Doppler parameter of 15kms _1 . Here, we are able to demonstrate that the 
power-law behavior persists for optically thin systems as well, in agreement with 
Keck measurements of the forest (Hu et al. 1995). 

The simulations reveal strong correlations between the cloud physical prop- 
erties and column density (Zhang et al. 1996a, b). As shown in the right panels 
of Figure 1, the higher column density systems tend to be dense and warm, 
becoming more rarefied and cool toward lower values. Most of the baryons in 
the universe, 60%, are found to reside in the optically thick absorbers. This 
fact alone is able to reconcile the BBN estimate for the density of the IGM with 
the low Gunn-Peterson limits and a QSO-dominated ionizing UV background, 
since the Lya opacity of a homogeneous component of the IGM would now be 
reduced by more than a factor of 5. 

The IGM, however, is found to fragment into even lower column density 
systems (Hu et al. 1995). While this may not be surprising at first glance, a 
comparison with the minimal column density that could result from Jeans insta- 
bility reveals a problem. A cloud with an internal density pb that is optically thin 
at the Lyman edge and in ionization equilibrium with the metagalactic radiation 
field will have a neutral hydrogen density of nni — 5.0 x 10 -15 cm -3 (pb / pb) 2 : (1 + 
z) 6 T 4 ~°' 75 r iI j _ 12 , where pb is the cosmic mean baryon density, here taken to 
correspond to fi&/i| = 0.05. The baryonic Jeans length in a medium of mixed 
dark matter and baryons is Xj = 27r(2/3) 1 / 2 c s (l + z)/H(z), where c s is the 
sound speed of the baryons associated with a linear perturbation and H(z) is 
the Hubble constant at redshift z (e.g., Bond & Szalay 1983). For an isothermal 

perturbation, Aj ~ 1 Mpc fy^T^l + z)~ 1/2 , or about 600 kpc at z = 3. This 
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gives a minimal column density due to Jeans instability of logAni — 13.5 at 
z = 3. The collapse of the system would result in an even larger neutral column 
density because of the increase in the recombination rate due to the increased 
density. Supposing that discrete features may arise only from Jeans unstable 
systems, one might then expect that lower column density systems cannot form. 
The detection of the low column density systems becomes even more perplexing 
on noticing the internal cloud densities required. Observations of the neigh- 
boring lines-of-sight to the images of the lensed QSO candidate HE 1104-1805 
(Smette et al. 1995), place a lower limit of 100 kpc on the size of the systems 
with logiVni > 13.2. For z > 3.5, such systems will be underdense, pb/pb < 1- 
The existence of the optically thin systems thus poses a puzzle. 

Figure 1 shows that in the simulation, optically thin features do form, and 
they indeed originate in clouds that are underdense. The fragmentation of the 
IGM is complete. At most a few percent of the baryons remain outside discrete 
systems to comprise a homogeneous diffuse component. Insight into the origin 
of the optically thin systems is provided by the simulation. I first discuss, how- 
ever, a simple model that offers a possible resolution to the difficulty of their 
occurrence. 

5. Cosmic Minivoids 
5.1. Spherical Model 

In an open universe, the growth of linear density perturbations "freezes out" at 
an early stage, at an epoch given roughly by 1 + zj « SI -1 — 1 (Peebles 1980). 
While bound systems will continue to collapse, unbound systems will expand 
like the expansion of the universe. In a flat universe, similar behavior will occur 
in a region that is locally underdense, since any perturbation inside will find 
itself in a locally open universe. 

Consider now a region that is slightly underdense in an Einstein-deSitter 
universe, a cosmic void. Because it is underdense, it will grow in size slightly 
faster than the surrounding cosmic expansion, increasing its underdensity. Ini- 
tially ripples within the void may grow, just as those outside, but once the void 
becomes sufficiently underdense, their growth will slow. Those that are bound 
will continue to collapse as they become nonlinear. But ripples that were orig- 
inally underdense, but not as strongly underdense as the background void, will 
freeze (relative to the cosmic background) into discrete structures within the 
void. These underdense systems would give rise to discrete optically thin ab- 
sorption features after exposure to an ionizing radiation field, even though they 
are unbound. 

This is illustrated in Figure 2. An initially underdense compensated spher- 
ical perturbation of comoving radius 3 Mpc is shown expanding in an Einstein- 
deSitter universe, with h^o = 1 and = 0.08. The evolution of the dark matter 
and the baryons is integrated as in Meiksin (1994). At the center of the under- 
dense region the density is initially raised slightly above the level of the void (in 
a growing mode), but still below the cosmic average density. A photoionizing 
radiation field is turned on at z = 6, increasing slowly with time (model ED of 
Meiksin & Madau 1993). A discrete absorption feature results along a line-of- 
sight passing through the void center, with a column density of logAni = 13.1 
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Figure 2. (Left panel) The evolution of density in a spherical 
minivoid. Upper curves show the dark matter density and the lower 
curves show the baryon density, both normalized by the cosmic mean 
density. The curves are shown at z = 3 (solid), z = 3.5 (dashed), 
and z = 4 (dot-dashed). (Right panel) The spectrum for a line-of-sight 
passing through the center of the minivoid. Optically thin absorption 
systems develop both at the minivoid center and at its boundary. The 
velocity is in the observer's frame. The solutions demonstrate that dis- 
crete absorption features may arise from a region that is Jeans stable. 

at z = 4. The perturbation is nearly frozen, with only a small decrease in the 
baryon overdensity due to pressure gradients, despite the fact that the system 
is Jeans stable. The reason is that the expansion velocity of the void from its 
center to the edge of the feature somewhat exceeds the isothermal sound speed 
of the gas. The feature is transient, but its expansion is dominated by the flow 
in the void rather than by thermal pressure. A second feature appears near the 
void boundary. This feature results from the collisional nature of the gas, which 
is compressed near the wall of the void where it meets the more slowly expanding 
external IGM. (This results in a positive bias for the baryons just within the 
void boundary.) Both mechanisms give rise to optically thin absorption features. 

5.2. Simulation Results 

The numerical simulation displays elements similar to this simple model. The 
underdense gas is indeed expanding faster than the cosmic rate, as shown by 
its positive peculiar velocity divergence in Figure 3. The structures form in 
regions a few comoving megaparsecs across. These minivoids are bordered by 
the collapsed filaments and sheets that give rise to the optically thick systems 
(Figure 3). By contrast, almost all of the He II systems, even those that are 
saturated, form in the minivoids. Because the effective opacity due to line- 
blanketing is dominated by features that are just becoming saturated (Madau & 
Meiksin 1994), the He II Lya opacity in a QSO spectrum serves as a probe of the 
structure of the minivoids. A more complete description of the He II absorption 
and the cloud structure is provided in Zhang et al. (1996a, b). 

These results admit a nontrivial consistency check based on the total opacity 
of the optically thin systems. While equation (1) may not be inverted to solve 
for Qd in the presence of clumping, a lower limit to the cosmic baryon density 
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Figure 3. Contour plots of a 150 kpc thick slice of the simulation box 
at z = 3, centered on a minivoid. The box size is 9.6 comoving Mpc. 
(Left panels) The baryon (top left) and dark matter (top right) contour 
levels are 0.5 (dotted), 1 (dashed), 3 (thin solid), and 5 (thick solid). The 
temperature contours, in units of 10 3 K, are 6 (dotted), 10 (dashed), 
14 (thin solid), and 20 (thick solid). The peculiar velocity divergence 
contours are, in units of the Hubble constant at z = 3, 5 (dotted), 
(dashed), -3 (thin solid), and -15 (thick solid). The low density region 
bordered by the filaments exhibits density fluctuations with positive 
velocity divergence. These fluctuations give rise to discrete optically 
thin absorption features. (Right panels) Contour plots of the HI and 
Hell column densities. The contour plots are at logiVHi = 12 (dotted), 
13 (dashed), 14 (thin solid), and 15 (thick solid), and at logA^Heii = 14 
(dotted), 15 (dashed), 16 (thin solid), and 17 (thick solid) (although 
radiative transfer has been ignored for Hell). The minivoid is populated 
by optically thin HI systems, and saturated Hell systems. 



may be obtained when the absorbers are optically thin and underdense. Taking 
T a to refer to the net stochastic absorption of the optically thin systems, £Id may 
be interpreted as the spatially averaged rms baryon density of the optically thin 
systems. (This is true only for optically thin clouds. See Jenkins & Ostriker 1991 
for the case including optically thick absorbers.) If the internal density of the 
clouds is less than the cosmic mean, provides a lower bound to the average 
cosmic baryon density. Counting the lines in Hu et al. (1995) for which the line- 
center Lya opacity is less than 0.5, and doubling to account for incompleteness 
(as discussed in Hu et al.), gives r Q ~ 0.08 for the optically thin sytems. This 

—3/2 

yields a lower limit for the cosmic baryon density of 0,f, > 0.02/i 50 , consistent 
with the value = 0.06 adopted in the simulation. 

The trends in Figure 1 permit a determination of the mass distribution of 
the clouds. Over the column density range 12.5 < logA^i < 14.5, the internal 

1/2 

baryon density of the clouds scales according to pb/pb — ^hii3> wnere AThi,13 
is the column density in units of 10 13 cm~ 2 . The neutral fraction scales like 
/hi = n m /n H ~ 3.8 x 10~ 6 [(1 + z) /AfT^ 7b Y^_ l2 NH^ lz . (Note that these 
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relations imply a roughly constant size for the absorbers of ~ 100 kpc over 
this column density range.) The observed number of absorption systems is 
d 2 N/dN m dz «6x 10 _13 JV^ 5 3 (1 + z) 2 - 5 (Hu et al. 1995). These results may 
be combined to derive the contribution of the Lya forest to the closure density, 

1 Amu H f, AT N m d 2 N 5/2 ( N m>msx \ 
tt Lya = / dN m-^ ~ AT a (1 + z r oclog — : . 2 

Pcrit C J /HI 3N m dz y iVffl.min / 

Most of the baryons lie in the column density range 12.5 < logA^ni < 14.5, 
distributed equally per decade in column density, in accordance with Figure 1. 
Only a small fraction, less than 5%, is found to reside in systems with column 
densities smaller than 10 12 - 5 cm~ 2 . Since the density of neutral hydrogen scales 
like the square of the baryon density, according to this model the reason it has 
been so difficult to detect the Gunn-Peterson effect is that the Lya opacity asso- 
ciated with a homogeneous component of the IGM has been suppressed by at least 
two orders of magnitude compared with the BBN value due to fragmentation. 

6. Observational Verification 

To date, most of the evidence supporting the results of the simulations is cir- 
cumstantial. The strongest evidence is provided by the large cloud sizes inferred 
from neighboring lines-of-sight to QSOs (Smette et al. 1995; Fang et al. 1996), 
however the filamentary nature of the absorbers found in the simulations has yet 
to be demonstrated. The possibility that the clouds are giant spheres (though 
perhaps unlikely), cannot yet be excluded. Indeed, the higher column density 
systems 10 16 cm -2 ) may still well arise from moderate sized spheroidal mini- 
halos (Rees 1986; Ikeuchi 1986). (See Charlton et al. 1996 for a discussion of the 
observational implications of the simulations for neighboring line-of-sight statis- 
tics.) A strategy for demonstrating a flattened or filamentary geometry for the 
absorbers is to map their distribution on the sky using multiple lines-of-sight 
with transverse separations up to a few megaparsecs. Similarly, while the sizes 
and abundances of the absorbers are consistent with their being the dominant 
reservoir of the baryons at high redshift (z > 2), the uncertainties in the cloud 
sizes and geometry, BBN baryon density, and magnitude of the ionizing UV 
background preclude a definite statement. 

Evidence in support of the findings of the simulations may come from an 
entirely different direction. Tytler et al. (1995) and Songaila & Cowie (1996) 
report the detection of C IV in the Lya forest. Songaila and Cowie find a 
median ratio of iVciv/^Hi ~ 3 x 10" 3 for 15 < logiV HI < 17. The detection 
of C II would permit a direct estimate of the ionization parameter U = n^/nn, 
where n 7 is the number density of ionizing photons and «h the total internal 
hydrogen density of a cloud. The Haardt & Madau (1996) spectrum has a 
large break at the He II Lyman edge, and a spectral index of ss 0.5 between 
the hydrogen and helium Lyman edges. Thus n 7 ~ 1.8 x 10~ 5 Fhi,-i2- Using 
the relation above found between cloud internal density and column density, 

—1/2 

the ionization parameter for the clouds is U — 2.3iV HI 13 at z = 3. There is 
some question as to how close the absorbers are to photoionization thermal 
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equilibrium. Low column density clouds will have too low a density (tih < 
10 _4 cm~ 3 ), to maintain photoionization thermal equilibrium, and so will be 
too cool (due to the adiabatic expansion of the IGM from which they arose), 
while infalling material in higher column density systems may be shock-heated 
to too high a temperature (Meiksin 1994). The temperatures are available from 
the simulation, and so may be used directly (Zhang et al. 1996b). Here we 
give indicative values for the metal ionization using CLOUDY (Ferland 1993), 
for which photoionization thermal equilibrium is assumed. Consider a typical 
cloud with a column density of logiVHi = 15.5. Its internal hydrogen density 
will be nn ~ 1 X 1CP 4 cm~ 3 and its ionization parameter will be U = 0.13. 
The neutral hydrogen fraction in the cloud is /hi ~ 10 -4 ' 8 , and the fractions 
of C IV, C II, and Si IV are / C rv ~ 10~ L4 , feu ~ 10~ 3 - 9 , and / S iiv ~ 10~ 52 . 
If the slabs are uniform, so that the metal column densities may be added for 
comparing with the H I, then matching to Aciv/Ahi ~ 0.003 requires a carbon 
abundance of about 0.4% solar. Both C II and Si IV will be undetectable, with 
Acn/Aciv ~ 3 x 10~ 3 and Asiiv/Aciv ~ 2 x 10~ 5 (for a solar abundance 
ratio of Si to C). A detection of C II matching the predicted value, however, 
may be made possible by stacking spectra showing C IV absorption. For a 
system with logA^i = 16.5, the situation improves. The simulations predict 
U j 0.04, which gives f m « lO" 43 , /civ ~ lO" ' 7 , feu ~ 10~ 21 , and / Si iv ~ 
10 2 ' 9 . The required carbon abundance is 0.2% solar. In this case, Nch/Nqw ~ 
0.04, so that C II should be detectable, while N S iw/N C Tv ~ 7 x 10~ 4 . It is 
critical that measurements be performed at high resolution to ensure that the 
same subcomponents of an individual Lya cloud are compared. Songaila <fe 
Cowie (1996) report values for Nqh/Nqiy in good agreement with the log iVm = 
16.5 prediction. Their values for Nguy/Ncw, however, are substantially higher 
than the estimates provided here. This suggests that there may be a second 
population of absorption systems not accounted for by the simulations (e.g., York 
et al. 1986), or that the assumption of photoionization thermal equilibrium is 
violated. A more careful calculation of the expected column density ratios using 
the simulation results must be performed to address this issue. 

While much remains to demonstrate the correctness of the structure of 
the IGM found in the numerical simulations, the ease with which the principal 
observed properties of the Lya forest are reproduced in a CDM-dominated cos- 
mology is a compelling argument in its favor. An additional prediction of the 
model emphasized here is that the Gunn-Peterson effect should be absent as 
a consequence of the fragmentation of the baryons and dark matter in cosmic 
minivoids. In the future, it may be possible to test rival cosmological theories 
based on the properties of the IGM. 

Acknowledgments. The author thanks his collaborators Peter Anninos, 
Mike Norman, and especially Yu Zhang, whose numerical simulation results 
provided the material for much of this talk, and to the conference organizers for 
their generous hospitality. 
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